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We report on two-particle azimuthal angle correlations between charged hadrons at for- 
ward/backward (deuteron/gold going direction) rapidity and charged hadrons at mid- rapidity in 
deuteron-gold (d + Au) and proton-proton (p + p) collisions at ^/snn = 200 GeV. Jet structures are 
observed in the correlations which we quantify in terms of the conditional yield and angular width of 
away-side partners. The kinematic region studied here samples partons in the gold nucleus carrying 
nucleon momentum fraction x ~ 0.1 to ~ 0.01. Within this range, we find no x dependence of the 
jet structure in d + Au collisions. 
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Observations in d + Au collisions at ^Jsnn = 200 GeV 
at the Relativistic Heavy Ion Collider (RHIC) reveal a 
significant suppression of hadron production at forward 
rapidity (deuteron-going direction) relative to p + p re- 
actions scaled up by the equivalent number of nucleon- 
nucleon collisions (N coU ) 0, 0- This suppression is 
observed for hadrons with momentum transverse to the 
beam direction over the range px ~ 1.5 — 4 GeV/c. In 
contrast, measurements at mid-rapidity 0, IE IE an d 
backward rapidity 0, 0, Q show a modest enhancement 
of the particle yield relative to N co u scaling over the 
same px range. Particle production at forward rapidity 
is sensitive to partons in the gold nucleus which carry a 
small nucleon momentum fraction (small Bjorken x).The 
suppression has generated significant theoretical interest 
including different calculational frameworks for under- 
standing the data [E H EE EH 

One such framework, the Color Glass Condensate 
(CGC), attempts to describe the data in terms of gluon 
saturation At small x the probability of emitting an 
extra gluon is large and the number of gluons grows in 
a limited transverse area. When the transverse density 
becomes large, partons start to overlap and gluon-gluon 
fusion processes start to dominate the parton evolution 
in the hadronic wave functions. As a result the gluon 
density becomes saturated. Since the nonlinear growth 
of the gluon density depends on the transverse size of the 
system, the effects of gluon saturation are expected to set 
in earlier (at higher Bjorken x) for heavy nuclei acceler- 
ated at ultra-relativistic energies than for free nucleons. 

In the leading order pQCD framework, a quark or 
gluon jet with large transverse momentum produced in a 
hard scattering process (high momentum transfer or large 
Q 2 ) must be momentum balanced by another quark or 
gluon jet in the opposite direction but with almost the 
same transverse momentum. Thus the azimuthal angle 
correlation between particles from the pair of jets (re- 
ferred to as di-jets) is characterized by two peak struc- 
tures separated by 180 degrees. In CGC calculations, the 
momentum to balance a jet may come from a large multi- 
plicity of gluons in the saturation regime, and thus no sin- 
gle partner jet may appear on the opposite side |l2|. This 
effect is analogous to the nuclear Mossbauer effect, and is 
often referred to as the appearance of mono-jets. Alter- 
native calculations, describing the suppression of single 
hadrons at forward rapidity in d + Au reactions in terms 
of leading twist pQCD effects, predict no such mono-jet 
feature (l3| . 

We want to probe this high gluon density regime in 
d + Au collisions with relatively high transverse momen- 
tum particles at forward rapidity. Such particles are 
likely to result from hard-scattering collisions involving 
small x partons in the gold nucleus. At small x the gluon 
density increases rapidly with Q 2 and saturation effects 
may be relevant for x ~ 0.01 at modest pt- CGC calcu- 
lations E3 predict significant suppression of the condi- 



tional yield and widening of away-side jet azimuthal cor- 
relations between rapidity-separated hadron pairs when 
one of those hadrons is at forward rapidity. 

In this Letter we report on measurements of two- 
particle azimuthal angle correlations between unidenti- 
fied charged hadrons in p + p and d + Au collisions at 
y/SNN — 200 GeV. In our analysis, the two particles are 
referred to as the trigger and associated particles. The 
trigger particle is at forward ( 1.4 < rj < 2.0) or backward 
(—2.0 < 77 < —1.4) rapidity and the associated particle is 
at mid-rapidity, |?y| < 0.35. The particles are separated 
by an average pseudorapidity gap < A?7>~ 1.5. The cri- 
teria for trigger particles, associated particles and event 
selection are described elsewhere |3j,[l4|. The two-particle 
azimuthal angle correlation technique has been used ex- 
tensively by RHIC experiments and is described in detail 
elsewhere [H EE EE EE EH ■ In this technique the az- 
imuthal correlation function is formed from the angular 
difference, A<p = (f> assoc — (f> trM J ; between each trigger and 
associated particle pair. Two jet peaks are normally ob- 
served in such correlation functions: the near-side peak 
(Acf> ~ 0) in which the two particles come from the same 
jet, and the away-side peak (A</> ~ it) in which they come 
from the back-to-back jets. In addition to these peaks 
the correlation functions usually also have a Aip inde- 
pendent combinatoric background contribution which is 
due to trigger-associated pairs from different jets or from 
non-jet processes. 

We can construct separate correlation functions that 
are sensitive to partons in the gold nucleus with differ- 
ent Bjorken x ranges. By choosing trigger particles with 
1.0 < pt < 5.0 GcV/c at forward (backward) rapidity, 
we sample partons in gold nuclei with x ~ 0.01(0.1). At 
x ~ 0.1, we do not expect to be sensitive to anysatura- 
tion effects, but we may be sensitive at x ~ 0.01 |l9|. The 
comparison in d + Au reactions between these two cases, 
as well as with the p + p case, may give insights into 
possible saturation effects on jet production and other 
mechanisms for forward rapidity single particle suppres- 
sion. It should be noted that the prediction of mono-jets 
in E2 assumes one particle at pseudorapidity 77 = 3.8 
and one at mid-rapidity, thus demonstrating sensitivity 
at even lower xau (~ 10~ 4 ) which is outside the range 
presented in this analysis. 

Data for this analysis were collected by the PHENIX 
experiment [2(J in the 2003 RHIC p + p and d + Au run- 
ning period. In the case of d+Au collisions, we divide the 
data into two centrality (impact parameter) classes based 
on the number of hits in the backward-rapidity PHENIX 
Beam-Beam Counter (BBC, -3.9 < 77 < -3.0). Cen- 
tral (peripheral) collisions comprise — 40% (40 — 88%) 
of the minimum bias cross section respectively. Using a 
Glauber model Q and a simulation of the BBC, we de- 
termine < N co u >= 4.7±0.4 for peripheral collisions and 
< N co u >= 13.0 ± 0.9 for central collisions. 

The trigger particles are measured in the PHENIX 
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muon spectrometers |20j. In this analysis we only se- 
lect trigger particles from 1.4 <| 77 |< 2.0 to obtain 
homogenous acceptance in transverse momentum from 
1< pr < 5 GeV/c and to reduce beam correlated back- 
grounds. We identify hadrons, as opposed to muons, 
in the muon spectrometers by comparing their momen- 
tum and penetration depth. This hadron identifica- 
tion method is described elsewhere 3]. It is notable 
that our trigger hadrons have a modified composition 
(pion/kaon/proton ratio) relative to that at the collision 
vertex due to species-dependent nuclear interaction cross 
sections. Detailed simulations show that kaons make up 
65 — 90% of positively charged trigger particles and pions 
make up 70 — 90% of negatively charged trigger particles. 
The baryon contribution to our trigger particle sample 
is negligible. We find the two-particle azimuthal angle 
correlations for positively and negatively charged trigger 
particles to be consistent and therefore combined the re- 
sults. The associated particles are unidentified charged 
hadrons measured in the PHENIX central spectrome- 
ters [2(| which cover | rj |< 0.35 and in this analysis 
have 0.5 < pr < 2.5 GeV/c. Standard track selection 
criteria 01 are applied. 

For comparison we have also included measurements 
where trigger particles and associated particles are both 
measured in the PHENIX central spectrometers at mid- 
rapidity. The d + Au points for this comparison are 
from | 14| and the p + p point is an extension in pt of 
the analysis that was published in [lij . 

We define the azimuthal angle correlation function as: 



2.0<Pt 9 <5.0 GeV/cO^pT^I GeV/c 



CF 



dN(A(/))/d(A(f ) ) 
acc(A(p) 



(1) 



where dN (A(f) / d(A(f) is the measured two-particle dis- 
tribution and acc(A(j)) is the two-particle acceptance ob- 
tained by an event mixing technique in which we mix 
trigger particles with associated particles from different 
events within the same centrality and collision vertex cat- 
egory. This correction is necessary because the PHENIX 
central arm detector is not azimuthally symmetric and 
the pair acceptance varies as a function of Ac/). 

In Figure n we show the correlation functions for trig- 
ger particles with pt = 2 — 5 GeV/c and associated par- 
ticles with pt = 0.5 — 1.0 GeV/c. The top panel is for 
p + p collisions where we have combined the results from 
forward and backward pseudorapidity since the collision 
system is symmetric. The middle and bottom panels of 
Figure^show the correlation functions for central d+Au 
collisions. The middle panel is for the trigger particle at 
forward rapidity and the bottom panel is for the trig- 
ger particle at backward rapidity. A clear peak is seen 
near A</> = tt in all cases corresponding to the away-side 
jet. It is notable that there is no peak near A<f) = 0, 
as expected, because the rapidity gap between the two 
particles is larger than the width of the near side jet. 
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FIG. 1: Azimuthal angle correlation functions. On the plots, 
the Gaussian widths from the fits and the signal to back- 
ground ratio integrated over n — 1 < A<j> < n + 1 are shown. 
Note that the t/-axis is zero-suppressed on the middle and 
bottom panels. 



After constructing the correlation functions in various 
bins in p T ssoc 1 p^ 19 and rf r%9 we used two methods to 
determine the unnormalized number of trigger-associated 
particle pairs, N pa i r , above a constant background. In 
the first method, we define 



ir+l 



N.„ 



J2 CF(A^)- J2 CF(A<P), (2) 



A(A=7T-1 



A0=-1 



where the first term is the integral of the correlation func- 
tion in the area of the correlation peak (tt — 1 < A</> < 
tt + 1) and the second term is the integral away from 
the peak (— 1 < Acp < 1). In the second method we 
fit the correlation function with a Gaussian distribution 
centered at A4> = tt plus a constant background. The 
values of N pa i r obtained by each method are found to 
be consistent and the small differences are included in 
our systematic errors. The solid lines in Figure ^ show 
the resulting fits and the Gaussian width parameters (a) 
together with the integrated signal to background ratios 
(Jj) over the signal region (tt — 1 < Acp < tt + 1) are 
quoted. 

The conditional yield (per trigger particle ) is defined 
to be 



CY 



N pair I £ as 



N, 



trig 



(3) 



where e assoc (~ 0.15 ± 0.015) is the efficiency times ac- 
ceptance for associated particles and N tr i g is the number 
of trigger particles used to generate the correlation func- 
tion. Eassoc is obtained for each colliding system, central- 
ity class, and transverse momentum bin by a GEANT 
based simulation of the PHENIX detector Q . 
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2.5<p" 9 <4.0 GeV/c, 1 .0<tf? sac <2.5 GeV/c 
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FIG. 2: Conditional yields are shown as a function of trig- 
ger particle pseudorapidity. The data points at mid-rapidity 
for d + Au collisions are from [T^| . To increase visibility, we 
artificially shift the data points belonging to the same rf rt9 
bin. The errors on each points are statistical errors. The 
black bar around 0.1 on the left of the plot indicates a 10% 
common systematic error for all the data points due to the 
determination of associated particle efficiency. There is an 
additional +0.037 systematic error on the mid-rapidity p + p 
point from jet yield extraction, which is shown as the arrow 
on that point (similar analysis as |T^|'). 



It is interesting to plot the conditional yields as a func- 
tion of r\ tT%9 . Changing rf rig from —2.0 to 2.0 effectively 
changes the range of Bj or ken x of sampled partons in gold 
nuclei from Q-ll lo4 to 0-0lio\oo7- Results are shown in 
Figure [21 The first observation is that there is no dif- 
ference beyond statistical fluctuations in the conditional 
yields for p + p, d + Au peripheral, or d + Au central 
collisions at any trigger particle pseudorapidity. 

We further quantify any nuclear modification in the 
conditional yield by defining the following ratio: 



IdAu 



CY 



d+Au 



CY 



(4) 



p+p 



Figure [21 shows the ratio IdAu as a function of p9j? soc for 
central and peripheral d + Au collisions, two different 
p^* 9 ranges, and for forward and backward rapidity trig- 
ger particles. In the plot, the shaded bands on each of the 
data points are the systematic errors which are the differ- 
ences in N pa i r obtained from the two methods described 
above and are independent from point-to-point. There is 
also a point-by-point correlated ~ 2% systematic uncer- 
tainty in the centrality dependence of s assoc determined 
by embedding Monte Carlo tracks into real events. The 
size of this uncertainty is comparable to the width of the 
IdAu = 1 line. 

Our measurement of the nuclear modification of the 
conditional yield indicates that IdAu with the trigger par- 
ticle at forward rapidity (sampling low-x partons in the 
gold nucleus) and backward rapidity (sampling high- a; 
partons in the gold nucleus) are both consistent with one. 
There may even be some evidence of slight enhancement 
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FIG. 3: I dA u vs. p^ S30C for different centrality, p% ig and r) trig 
bins. To increase visibility, we artificially shift the data points 
belonging to the same py ssoc bin. 
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FIG. 4: The ratio of correlation peak widths between d + Au 
and p + p collisions. Only statistic error is shown. To increase 
visibility, we artificially shift the data points belonging to the 
same p?? soc bin. 



for the case with trigger particles at forward pseudorapid- 
ity in central d + Au collisions. We note that if mono-jets 
were a major contributor to our trigger particle sample 
in our x range, we would have expected a decrease in 
the conditional yield for d + Au central collisions when 
the trigger particle is at forward pseudorapidity. Our 
measurement is inconsistent with any large nuclear sup- 
pression (i.e. mono-jets) of the jet structure in this kine- 
matic range. However, we note that in these modest 
Pt ranges, there may be contributions from both "hard" 
(large Q 2 ) processes and "soft" coherent (small Q 2 ) pro- 
cesses. In d + Au collisions "soft" particle production 
is shifted away from forward rapidity towards backward 
rapidity [2J. Thus, the fraction of hadrons at forward 
rapidity from "hard" processes may be increased in cen- 
tral d + Au reactions, which may offer an explanation for 
the modest enhancement seen in the conditional yield for 
this case and could also mask off small mono-jet signals. 
We have also compared the Gaussian widths of the 
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correlation peaks in d + Au collisions to the widths in 
p + p collisions. The ratios of the d + Au widths to the 
p+p widths are plotted in Pigure|Uas a function of p^ ssoc . 
There may be a hint of a slight p9j? soc dependence in the 
ratio, but overall there is no significant difference in the 
width in d + Au collisions for different rf r%9 . 

In conclusion, we have measured the two-particle az- 
imuthal angle correlations in two centrality categories of 
d + Au collisions and in p + p collisions. We measured 
trigger particles at forward, backward and mid-rapidity 
and correlated them with associated particles at mid- 
rapidity. The associated particle conditional yields in 
central d + Au collisions are consistent with the condi- 
tional yield in p + p collisions. These conditional yields 
also do not change as we vary the trigger particle pseudo- 
rapidity over the range | rj |< 2.0. We have also compared 
the widths of the away-side jet peaks in d+Au and in p+p 
collisions, and find no evidence for large 77*" 9 -dependent 
modification. 
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